Possible DD and BB Molecular states in a chiral quark model 
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We perform a systematic study of the bound state problem of DD and BB systems by using 
effective interaction in our chiral quark model. Our results show that both the interactions of DD 
and BB states are attractive, which consequently result in I G (J PC ) = + (0 ++ ) DD and BB bound 
states. 
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(N ' I. INTRODUCTION 

t-H ■ 

o : 

\ Since the discovery of X(3872), many X, Y, and Z exotic states have been reported. These hadrons have ever been 
! * . explained as molecules, tetraquarks, hybrids et al. because they can't be interpreted as simple quarkoniums. In our 
previous chiral quark model calculation it is found that the newly observed hadrons, such as .2f&(10610), Zb(10650), 
O ! X(3872) and Y(3940), might be assigned as the BB* , B*B* , DD* and D*D* bound states. These results stimulate 
' our further interest in studying their analogues, i.e. the systems of DD and BB. 

So far, several works have been done to calculate the DD and BB states (2T4ll|. Valcarce et al. considered 
the DD coupled to charmonium-light two- meson systems, like J/ipuj channel, and concluded that the + (0 ++ ) DD is 
r— i the only possible bound state. Ke et al. @ supported the existence of + (0 ++ ) DD and BB molecules in the Bethe- 
'■^ . Salpeter framework. In one-meson-exchange model, Liu et al. calculated the binding energies of DD and BB 

systems. Meanwhile, Zhang et al. [1] got such molecules: BB of 10580±100 MeV and DD of 3760±100 MeV in QCD 
sum rule calculation on the quark level. Their results are consistent with Wong's prediction in a two-gluon-exchangc 
model [lo| . On the contrary, Yang et al. [ll[ argued that in color-singlet channel DD and BB molecules didn't exist. 
To sum up the above calculations, one sees that the existence and properties of the possible DD and BB molecular 
states are presently model dependent. Further theoretical investigations are expected to be significant. 
£NJ | In this paper, wc perform a dynamical study of the DD and BB systems with isospin 1 = and 1 in our chiral 
^ • quark model by using the effective interaction. The chiral quark model was built in such a way that the chiral 
symmetry is restored by introducing the coupling between quark field and Goldstone bosons and the constituent 
quark mass is obtained as a consequence of the spontaneous vacuum symmetry breaking. In our chiral SU(3) quark 
model, onc-gluon-exchangc (OGE) governs the short range and scalar chiral field as well as pseudoscalar chiral field arc 
induced for restoring the chiral symmetry. As is well known, the short range mechanism of the quark-quark interaction 
mechanism is still an open problem, and people is debating whether the OGE plays a dominating role in the short 



| range, or vector meson exchange docs, or both of them arc important. Thus to examine the short range mechanism, 
I ■ we developed our chiral SU(3) quark model into the extended chiral SU(3) quark model in which the vector meson 
• • ' exchange is included. During the past few years, both the chiral SU(3) quark model and the extended chiral SU(3) 
. — ' ci uar k model have appeared to be quite successful in reproducing the spectra of the baryon ground states, the binding 
k> ' energy of deuteron, the nucleon-nucleon (NN), kaon-nucleon (KN) scattering phase shifts, and the hyperon-nucleon 
^ ! (YN) cross sections [H-[ll|. 

Recently, the chiral quark model was also employed to study the interactions and structures of the heavy-quark 
systems by using the Resonating Group Method (RGM) [17H19| . Here we will use the same chiral quark model to 
study the DD and BB interactions. Different from the RGM method, we have derived analytical forms of the total 
interaction potentials between the two S-wave heavy mesons as discussed in Refs. [H, We thoroughly investigate 
the possible bound states of DD and BB systems by solving the Schrodinger equation with our analytical potentials 
between the two clusters, anticipating that this method would give a more accurate description of the short-range 
interaction between the two clusters than the RGM docs. 

The paper is organized as follows. In section [Til the framework of our chiral quark model is briefly introduced, 
and the analytical forms of the effective interaction potentials between the two S-wave heavy mesons under our chiral 
quark model are given. The bound state solutions for the DD and BB systems are shown and discussed in Sec. IIIII 
Finally, a short summary is given in Sec. IIVI 
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II. FORMULATION 



In our chiral SU(3) quark model, we consider the scalar meson-exchange interaction and the pseudoscalar meson- 
exchange interaction which arc induced from quark-chiral field interaction, OGE interaction and confinement potential. 
Moreover, in our extended chiral SU(3) quark model, we also take the vector meson-exchan ge in teractions into account 
which would almost replace the OGE interaction in the original chiral SU(3) quark model In principle, the axial 
vector meson exchange should also be included to keep the chiral symmetry. However, as the axial mesons are much 
heavier than the chiral symmetry breaking scale, they arc simply omitted in our calculation. 

The framework of our models has been discussed extensively in the literature (l2l - fl9l [2l| - |2^ . In this work, different 
from the RGM calculation, the internal kinetic energies and the internal interactions of each meson are not necessary 
to be calculated. As a result, the Hamiltonian of relative motion in this work reads 

H = T rel +V eff , (1) 

where T re i is the kinetic energy operator of the relative motion between the two mesons, and V e ff is the effective 
interaction potential derived from the quark-quark (quark-antiquark) interaction between two mesons by integrating 
the internal coordinates £1 and £2 of two mesons: 

v -ff = E / ri&)v*2&)v(?i i )vi(&)M&)<&&, (2) 

ij J 

while <pi(£i) and (^2(^2) are the intrinsic wavefunctions of two S-wave mesons, taken as one-Gaussian form: 

P(0 = (— ) e 2 ■ ( 3 ) 

Here, jj, is the reduced mass of the two quarks inside each meson and uj is the harmonic-oscillator frequency of the 
meson intrinsic wavefunction. V(fij) in Eq. 2 represents the interactions between the i-th light quark or antiquark 
in the first meson and the j-th light quark or antiquark in another. 

In this work, because there is no color-interrelated interaction between the two color-singlct clusters, such as OGE 
interaction and confinement potential, we only consider the meson-exchange interactions. Therefore, for the chiral 
SU(3) quark model 

Vifa) = E^"(^) + E^"(^), (4) 

a=0 a=0 

and for the extended chiral SU(3) quark model 

vifa) = E yaa fo) + E v * a fa) + E VPa (*u ( 5 ) 

a— a— a— 

with V aa (fij ) and V* a (fjj ) being the interactions respectively induced from scalar meson exchange and pseudoscalar 
meson exchange. V Pa (fij) indicates the vector meson exchange interaction. For quark- q uar k ( ant iq uark- ant iq uark) 
interaction, V' Ja (rij), V^ a (fij) and V Pa (fij) have been described in detail in Refs. [l|, [UtH, Erllll : 

V aa {n-j) = -C(g ch ,m aa ,A)X 1 {m aa ,A,r ij )(\?\'*), (6) 

O 

Tfl 

V^inj) = C(g ch ,m„ a ,A) ff ° X 2 (m 7r , A, r^-) x • a,) (A"A°) , (7) 

lzniimj 



V Pa {nj) = C(g chv ,m Pa ,A) 



Mrn Pa ,A,r ij) + ^x(l + ^^^ + ^^ 
Qmitrij V gchv M N g 2 chv M N 



x X 2 (m Pa , A,r 4j ) fa ■ Oj) 



(A?A|) , (8) 



with 



,2 A 2 



C(g ch ,m,A) = ~ , m, (9) 
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Xi(m,A,rij) 
X2(m,A,nj) 
Y(x) 



Y(m nj ) - ^Y(A nj ), 
Y(mr t3 )-f±) Y(Ar l3 ), 



— e 

x 



(10) 

(11) 
(12) 



where A a is the Gell-Mann matrix in flavor space, and A is the cutoff mass which indicates the chiral symmetry 
breaking scale, uii and rrij are the masses of the i-th light quark or antiquark in the first meson and the j-th light 
quark or antiquark in another respectively, while m aa , rn^ a and m Pa in Eqs. [6j [7j and [8] are the masses of the scalar 
nonets, the pseudoscalar nonets and the vector nonets, respectively. Mn in Eq. [5] is a mass scale usually taken as 
the mass of nucleon [l4j]- g c h is the coupling constants for the scalar and pseudoscalar nonets. g c h v and f c h. v are the 
coupling constants for the vector coupling and tensor coupling of vector nonets, respectively. 

In DD and BB systems, we don't consider the one-meson exchange interactions between two heavy quarks or 
between one heavy quark and one light quark, because if one wants to include the interactions related to heavy 
quarks, the heavy meson-exchanges as well as light meson exchanges must be also considered simultaneously, and 
these interactions arc beyond our SU(3) models. By using the method described in Refs. 0, Ho[ and integrating the 
internal coordinates of two mesons, we get the analytical effective interaction potentials between two mesons DD(BB) 
as 

8 8 8 

V efJ (R) = £ V£ (R) + £ Ki (£) + E Ki (R)> 



with 



V q f(R) 
V q P i(R) 



-Ga a C{g ch ,m aa , A)X lq g(m aa , A, R) (A£A£) , 

o 

Tfi 

G^ a C(g ch ,m WcL ,A) * a X 2 qq(m na ,A, R) x (a q ■ a q ) {\ q X q ) 



(13) 
(14) 



G Pa C(g chv ,m Pa ,A) 



X lqq\ m Pa , A ,R) + 7 1 H 77 r — 772- 



xX 2qq (m Pa , A, R) (<r q ■ a q ) 



( ■ 



(15) 



Here, G aa ^ a . Pa is the G-parity of the exchanged meson, and R is the relative coordinate between two different mesons, 
namely, the relative coordinate between the two centers-of-mass coordinates of the two mesons, and 



X lqq (m,A,R) = Y qq (mR) - -Y qq (AR), 

m 

X 2qq (m,A,R) = Y qq (mR)-(-) Y qq (AR) 



(16) 
(17) 



In above equations, m q and m q are masses of the light quark and antiquark, respectively. The modified Yukawa term 
in Eqs. [TBI and [T71 reads 



Y qq (mR) 



Here, 



I mq \ ' 
\ rriQ+niq J 



m q +m( 



(18) 



(19) 



rriQ and rriQ are masses of the heavy quark and antiquark respectively, and fi q Q = m " +r ^ Q ■ 

There are some necessary parameters in the potentials of our chiral quark model. In this work, we adopt the 



parameters determined in our previous works [2l|-[23|. The up/down quark mass m q is fitted as the nucleon 
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mass and taken as Mjv/3 ~ 313 MeV. The coupling constant for the scalar and pseudoscalar chiral fields g c h = 2.621 
is fixed by the relation of 

9ch _ 9 9nNtt m u 



4tt 25 4tt M%' 

with g^ N7r /4ir = 13.67 determined from experiments. 

In our extended chiral SU(3) quark model, the vector coupling constant g c h v and tensor coupling constant f c h v hi 
Eqs. [U [15] are fitted by the mass difference between TV and A, when the strength of the OGE is taken to be almost 
zero. When the tensor coupling is neglected, g c hv = 2.351 and f c hv = 0; when the tensor coupling is considered, 
g c hv = 1-973 and f c h v = 1-315. The harmonic-oscillator frequency w, equal to l/(m„6„) = l/{m c b1) where b u is fitted 
by the N — N scattering phase shifts, is taken as 2.522/m" 1 in the chiral SU(3) quark model and 3.113/m 



the extended chiral SU(3) quark model. In our calculation, the masses of the mesons are taken from the PDG [25jj, 
except the a meson, which does not have a well-defined value. Here m a is obtained by fitting the binding energy of 
the deuteron [lij]. It is m G = 595 MeV in our chiral SU(3) quark model, 535 MeV for neglecting tensor coupling and 
547 MeV for considering tensor coupling in our extended chiral SU(3) quark model. The cutoff mass A is the chiral 
symmetry breaking scale and taken as 1100 MeV as a convention. 

The remaining parameters to be determined are the heavy quark masses m c and m&. In our work, we find that 
the final results are not sensitive to the variation of the heavy quark masses, and we take m c = 1430 MeV [26[ and 



m& = 4720 MeV [27| as typical values. 



III. NUMERICAL SOLUTIONS 



A. DD 



Here we study the DD system with different isospin /. Following the approach introduced in section II, we get the 



analytical effective interactions between DD in our chiral quark model which are depicted in Fig. 1(a) From Fig 



1(a) one sees that the DD interaction is attractive, especially for isospin 1 = case. To study if such an attraction 



is strong enough to bind the DD system, we solve the Schrodinger equation with the programs developed in Refs. 

[29| . and list the obtained binding energies in Table fl] Because the mass difference between DD and D*D* is 
not very large, we also consider the coupled effect of D*D* on DD. We carry out a perturbative calculation to see 
the contributions of the off-diagonal elements of this coupled channel to the binding energies of DD and list them in 
Table HI From Table QJ we see there is only one I G {J PC ) = + (0 ++ ) S- wave DD bound state with a binding energy 
3-35 MeV in our chiral quark model. 



TABLE I: The binding energy and the root of mean square radius of DD binding system. The binding energy is listed in such 
a way: B+AB, B is the binding energy deduced in the single channel calculation and AB is the perturbation correction value 
deduced from the off-diagonal elements of the coupled channel of DD and D* D* . 



X-SU(3) QM 


Ex. x-SU(3) QM 




gchv — 2.351, j c hv — 


g c hv = 1.973,/chi, = 1.315 


B+AB (MeV) r rms (fm) 


B+AB(MeV) r rms (fm) 


B+AB(MeV) w(fm) 


I = 1.0+2.3 3.7 


33.3+1.4 1.0 


21.6+0.2 1.1 


I = 1 







Further analysis shows that the DD interaction is dominated by er, cr', uj and p exchanges. For the 1=0 case, in 
the chiral SU(3) quark model, both a and a' exchanges provide attractive interactions, so the total interaction is 
strong enough to form a DD bound state. In the extended chiral SU(3) quark model, the contributions of vector 
meson exchange are also included, and p and u exchanges provide additional attraction. Therefore, the DD system 
has a larger binding energy as shown in Table |H We see that the perturbative contribution of the coupled channel is 
not very large compared to the results of the single channel and wouldn't obviously change the main feature of the 
binding solution. The radii of D and D both are about 0.54/to in the chiral SU(3) quark model and 0.49/m in the 
extended chiral SU(3) quark model (since uj is different in the two models), and we find that the rms radius of this 
+ (0 ++ ) DD bound state is bigger than the sum of the radii of D and D, so we conclude that this DD could form a 
molecule and this + (0 ++ ) DD molecular state has a mass 3695-3726 MeV. For the 1=1 case, in the chiral SU(3) 
quark model, a exchange provides attraction but a' exchange provides repulsion, thus the total attractive interaction 
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(a) 



(b) 



FIG. 1: The total interaction potential of DD and BB systems. 1(a) the total interaction of DD system and 1(b) the total 
interaction of BB system. The solid, dashed and dotted lines represent the results obtained from the chiral SU(3) quark model, 
the extended chiral SU(3) quark model without and with the tensor coupling of vector field, respectively. 



is too weak to make a DD bound state. In the extended chiral SU(3) quark model, the additional attraction provided 
by bj exchange and the additional repulsion provided by p exchange almost cancel each other, as a result, the total 
interaction of all the meson exchanges is still too weak to bind DD. Thus, no 1~(0 ++ ) DD bound state can be 
obtained due to the insufficiency of the attraction. 



B. BB 



TABLE II: The binding energy and the root of mean square radius of BB binding system. The binding energy is listed in such 
a way: B+AB, B is the binding energy deduced in the single channel calculation and AB is the perturbation correction value 
deduced from the off-diagonal elements of the coupled channel of BB and B*B* . 



X-SU(3) QM 


Ex. x-SU(3) QM 




g chv = 2.351,/ c(ll) = g chv = 1.973,/ ch „ = 1.315 


B+AB(MeV) r rms (fm) 


B+AB(MeV) r rms (fm) B+AB(MeV) r rms (fm) 


7 = 19.2+1.2 0.9 


81.8+11.3 0.6 64+2.6 0.6 


I = 1 





Considering the resemblance between DD and BB, the one meson exchange interaction potentials should have 
similar properties. In our calculation we find the total interaction potential of BB system is also attractive, as shown 



in Fig. 1(b) By Solving the Schrodinger equation and considering the perturbative correction of the coupled channel 
of B*B* to the binding energies of BB, we list our results in Table |TIJ We notice our results are similar to the ones 
of Liu et al. [l7| ■ It should be emphasized that the BB system is easier to form a bound state than the DD system, 
because BB system has a much heavier reduced mass 2640 MeV than that of DD system of 932 MeV. We also find 
that the radii of B and B are both about 0.52/m in the chiral SU(3) quark model and 0.46/to in the extended chiral 
SU(3) quark model. However, according to Table |TT] we see that the rms radius of this BB bound state is smaller 
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than the sum of the radii of B and B. This feature means B and B are overlapped with each other, and this + (0 ++ ) 
BB bound state seems not a molecular state but might be a tetra-quark state with a mass 10467-10540 MeV in our 
chiral quark models. 

IV. SUMMARY 

In this work, we have studied the bound state problem of DD and BB systems with an analytical effective interaction 
potential in our chiral quark model. We find only one + (0 ++ ) DD molecule could exist with a mass 3695-3726 MeV, 
agreeing with the predictions of Rcfs. 0413 G3- For the + (0 ++ ) BB system, they are bound even deeper than the 
DD, however, the small rms radii mean the overlap is sizeable and, therefore, it may not be interpreted as molecular 
states, but a tetra-quark state. 

We expect future experimental measurements would test the results of our chiral SU(3) quark model and the 
extended chiral SU(3) quark model. Since our calculation shows the important role of the vector meson exchange, we 
hope that the future measurement could recognize whether vector meson-exchange or one-gluon-exchangc dominates 
the short range interaction between quarks. 
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